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Abstract 

We study the change of entanglement under general linear transformation of modes in a bosonic 
system and determine the conditions under which entanglement can be generated under such trans- 
formations. As an example we consider the thermal entanglement between the vibrational modes 
of two coupled oscillators and determine the temperature above which quantum correlations are 
destroyed by thermal fluctuations. 
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1 Introduction 



Consider two ions in an ion trap, or two atoms in a solid vibrating around their equilibrium positions. 
We ask how much pure quantum correlation or entanglement exists between the vibrational modes 
of these two ions or atoms at a given temperature? How much we can raise the temperature before 
the quantum correlation between the ions is destroyed? In a more general context, we can ask the 
degree of entanglement of two bosonic modes in a many body system. This property known as thermal 
entanglement has been intensively studied mostly for the spin degrees of freedom of atoms or ions in 
the past two or three years ^[21 El- The reason for this restriction has been threefold. The first is that 
a calculable measure of entanglement when the two systems are in a mixed state has been known only 
for two dimensional systems jlj. The second reason has been obviously the interest in two dimensional 
systems as representatives of qubits in quantum computers 0. The third is that spin systems are the 
prototype of interacting many body localized fermions and it has been desirable to see what happens 
to entanglement when a spin system undergoes quantum phase transition 6 . 

While the pursuit of this problem for intermediate dimensions has been impossible (for the lack of 
a measure of entanglement) or of little interest, all the above three motivations exist also for the 
other extreme of dimensionality, namely systems with continuous degrees of freedom. First, for a 
class of states called symmetric Gaussian states, a closed formula exist for their entanglement of 
formation second, systems of continuous degrees of freedom are of wide interest as candidates for 
the implementation of quantum computing jHj and third, in condensed matter systems, the continuous 
degrees of freedom like the vibrational modes of localized atoms in a crystal or the bosonic modes of 
a system of identical particles can also be entangled. The same phenomenon is of interest in bosonic 
field theories [T3] . 

When dealing with systems of identical particles, it is the modes and not the particles which should 
be treated as subsystems [S] and one can quantify the entanglement of these subsystems in a second 
quantized approach [Hj IT U I12 | ITHj. In this case entanglement changes by redefining the modes, a 
non-local operation on subsystems. 

In this paper we study a system composed of identical spinless bosons subject to a free quadratic 
hamiltonian and study the thermal entanglement properties of different arbitrary modes of this system. 
This problem will also be of relevance when we want to understand how thermal fluctuations will affect 
the efficiency of protocols based on gaussian states [Ti] . 
So we consider a hamiltonian of the form 
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with [b a , bg] = 5 a /3 and a general linear mode transformation 
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To ensure the correct commutation relations between the new modes, the matrices S and T should 
satisfy: 

[S\T}=0, S ] S-TT^ = I. (3) 
We ask how much entanglement exists between modes k and I. This requires that we calculate the 
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reduced density matrix of these two modes, pti ■= tr^pth) where pth = 
k, I means that the modes k and I are excluded when taking the trace. 

We will show that the reduced density matrix of any two modes, is always a gaussian state(?). We 
obtain the condition under which the transformations (J2J) can not produce entanglement. We then 
consider an example which is a prototype of a wider class and calculate exactly the entanglement 
between the two modes and its dependence on temperature, specifically we obtain the threshold 
temperature above which quantum correlations are destroyed. 

We begin by collecting the necessary ingredients about gaussian states that we need in the sequel. In 
the Hilbert space of two harmonic oscillators a density matrix p is called a two mode gaussian state 
if its characteristic function, defined as 
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is a gaussian function. This can be written compactly as 

C{z* ll zi,zl,z 2 ) = e~ 
where zt = ( z* z\ z| z 2 ) an d M, parameterized as 
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where v 



is called the covariance matrix |15j . 

This matrix encodes all the correlations in the form M rs 
(a\, a±, a 2 , 0,2) and v := (01, a\, a 2 , a^) 7 . 

The conditions of separability of the two modes have been studied in a number of works |161 117j . In 
particular it has been shown that by a canonical transformation the covariance matrix can be put that 
when the covariance for a class of in which m c = mi = ttt-2 = 0. The conditions of separability then 
simplify to the following inequalities : 
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The above conditions only determine the separability of a gaussian state and not the amount of its 
entanglement. For a class of gaussian states invariant under the interchange of the two modes and 
called symmetric states, one can actually calculate in closed form the entanglement of formation 0. 
The closed formula is given in terms of the covariance matrix T defined as T rs = {r] r r] s + r) s r) r ), where 



r] := (x 1 ,pi,x 2 ,P2) , x r := -^(a r + a\) and p r := ^(4 " 
given in the covariance matrix T will have the form: 
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For symmetric gaussian states (those for which det a = det (3) one can apply local symplectic transfor- 
mation |181 119j and without changing their entanglement put their covariance matrix into the normal 
form 



(8) 



where k x > k p > 0. These new parameters are derivable from the symplectic invariants of the matrix 
M: 

n 2 = 4deta, k x k p = 4| det 7|, 
{n 2 - kl){n 2 - k 2 p ) = 4detM. 
The entanglement of formation of these symmetric states is then given by 

E f (p) := c+logc+ - c_ logc_, 

where c± = , and 

A = min ^1, 5 := \J {n — k x )(n — k p )^j . 
The state is entangled only when A < 1. 

Finally we note that for a state like p a := L e -^ab a b a w j iere fo a anc j are ^he usual harmonic 
oscillator operators, = 1, and Z = (1 — e~^ a )~ l , the single mode characteristic function 

C a (z*,z) := tr(e zb ^- z * b Pa ) is found to be 
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One way for doing this calculation is to expand the trace in the eigenstates of b^b. Then by using the 
properties of coherent states \z) := e zfet |0), (i.e. b\z) = z\z), (z\j) = e 2 * 7 and ^ J <i7d7*e~' 7 ' 2 |7)(7| = 
I), one writes 



(n\e- z * b e zb] \n) = —l-z\b n tf n \z) 
n! 
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By summing over n and performing the resulting gaussian integration one arrives at the stated result. 
Also for any two commuting modes 



C a f3{z\, Zl, Z2, z 2 ) — C a [Z 1 , Z\)Cp (z 2 , Z2). 
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This completes our short review of gaussian states. 

We now consider our mode transformations (|2"|). To express the covariance matrix of any two modes, 
say the modes k and I it is useful to introduce a compact notation. Let 
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and define the positive inner product between any two such vectors as 



(X\Y):=J2™th^X* a Y a . (15) 

a=l Z 

The characteristic function of the two modes is given by 

C kl (z* l ,z 1 ,z*,z 2 ) = tr (e z A-^+^l-^i ®L =1 Pa) , ( 16 ) 

By noting from (Q) that = p a , and inserting © in (J16|). rearranging the terms in the exponential 
and using (|11|) and (|13|) we find 



C M = e~5^=i \^S ka -z*T ka +z 2 S la ~z*T la \ 2 cosh ^ _ 

Comparing with (0J \^ we read the following matrix elements of the covariance matrix: (r = k, I): 

n r = -((S r \S r ) + (T r \T r )) , m r = -(S r \T r ), 
m s = ±((S k \Si) + (T^Tt)), 

m c = ^((SkM) + (S t \T k )). (17) 

We now prove that in any transformation which leaves the vacuum invariant, the new modes are 
disentangled. Any such transformation is one in which Tj a = Va. To show this we note that in this 
case the covariance matrix defined in © will have the parameters mi = m 2 = m c = 

111 

ni = -(S k \S k ), n 2 = -(Si\Si), m s = -(S k \S{) . 

Therefore in view of (j^J) the state will be separable if and only if 

(S k \S k ) > 1 , ((S k \S k ) - ms^Si) - 1) > |(S fe |^}| 2 . 

Due to (jHJ) we know that J2a=i \S ka \ 2 = 1 and J2a=i S ka Si a = 0. Since coth > 1, it is obvious that 
the first inequality is satisfied. If we know introduce a new inner product as 

(X\Y) :=f>*Y Q (coth^-l), 

Q=l L 

the second condition takes the form 

{S k \S k )(Si\Si) > \iSk\S0\ 2 , (18) 
which is satisfied by the Cuachy Schawrz inequality. This completes the proof. 

The above argument is also true for the transformations for which S = and T / 0. These new 
modes can also be called vacuum preserving by renaming the new creation and annihilation operators. 
The only transformations which can produce entanglement are those for which neither S nor T van- 
ish. Thus for systems of identical particles, these kinds of transformations play the role of non-local 
transformation which can produce entanglement. 
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2 An Example 



We will now consider such a transformation and for this purpose we choose an example which shows 
that our formalism is also applicable to systems of identical but localized and distinguishable bosonic 
particles. Consider two particles (atoms) oscillating around their equilibrium positions modeled by a 
mass-spring system with a Hamiltonian 

tt Pi , P2 , 1 2 , 1 2 1 1 2/ \2 

H = y + y + + + T^oOi - afc) , 

where x^ and denote the canonical coordinate and momentum of the i-th particle. Note that the 
particles are localized and distinguishable in this case. We want to calculate the thermal entanglement 
of these two atoms and study its dependence on temperature. In particular we want to see if there is 
a threshold temperature above which entanglement vanishes. This type of study has been intensively 
carried out for spins systems ^ |21 El HHj an d to our knowledge this is the first time where thermal 
entanglement of continuous degrees of freedom is being studied. 
We first find the the normal coordinates of the system; 

1 1 

Xt = — / =(x 1 +x 2 ), X 2 =: -j={xi- X2), 

p i ■= -/=0i +P2), P2 =■ -/=(pi -Pa), 



and the oscillator modes b\ := + iP\) and b 2 := '^ 2 ^( LJ X 2 + iPz), where u := y 1 + 2oJq. 

These modes diagonalize the Hamiltonian 

H = b\bi +ujb\b 2 , 

where we have ignored an overall constant. 

For calculating the thermal density matrix of the two particles, we proceed as before to determine 
the characteristic functions of the two modes where mode now means the degree of freedom of each 
individual atom. That is we have a\ = 772(^1 + Wi) an d a 2 = ~^\ x 2 + W2)- The relation of the 
new modes with the old ones turns out to be 



-L(b 1+li+ b 2 +^bi 



a 1 

02 = -j= [bi ~ £+b 2 ~ i~b 2 j , 
where £± = \ (^foj' 1 db y/u^j ■ From this equation and (|2I14|) we find that 

and 
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In order to simplify the notation lets us set u := coth f and v := coth ^p, a := i(o> 1 + and 



6 := i(w 1 — cj)u. With these conventions we will find from 



M = I 
4 
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Using the relations Q) we find the following form for the matrix T : 



u — a 



and jn|| that: 
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u — a + b 
u+a+b 

u+a-b j 

, from which we can calculate the 



This is not yet the final symmetric form of the matrix V as in 
entanglement. For this last step we need the parameters n, k x and k p which can be derived from the 
symplectic invariants of T. Actually it is simpler to do a canonical transformation Xi — ► axi, pi — > 
with a := (i^j^rf) 3 to put this matrix in the symmetric form (JHJ and read the parameters n, fc 2 and 



The result is: 
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From this last result and the definition of A in we find the condition of entanglement of the two 
atoms 

r2 



S z = u (a + b) 



coth | coth &± 



< 1. 
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Since u > 1, this inequality can be satisfied below a threshold temperature T c obtained by setting 
S 2 (T C , uj) = 1. Inserting the value of <5 from (|22|) in (|l(jp and then using Q we obtain the entanglement 
between the two atoms as a function of temperature and frequency. The result is plotted in figure 
1. The entanglement at zero temperature is obtained by inserting the value of delta in this limit, 
A = A=, in @. This leads to 



E(T = 0) = xlnx - (x - l)ln(x - 1) 



where x := ^77=^ ■ This is an increasing function of u as shown in figure 1. 

The maximum entanglement behaves like ( 161 2 + 1 — | log 2 (w — — l) 2 for small frequencies 
u ~ l(o;o << 1) and like ^ — 2 + | log 2 10 for large frequencies cj >> 1. 
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Figure 1: Entanglement of the two oscillators as a function of the natural frequency and temperature. 
The threshold temperature increases almost linearly with frequency. 

In summary we have set up an easy formalism for calculating the thermal entanglement of arbitrary 
bosonic modes in a rather large class of problems. For example the formalism can be applied to 
arbitrary lattices of bosonic modes to see how thermal fluctuations affect frustration of entanglement 
|21j . or to a chain of coupled oscillators. In the later case the entanglement can be obtained as a 
function of both the temperature and the distance between the particles. This is in contrast to the 
spin systems where due to the complicated nature of their spectrum this later dependence can not be 
obtained except only at low temperatures and under certain assumptions |2*2*] . 

I would like to thank V. Karimipour, A. Bayat, I. Marvian, L. Memarzadeh and A. Sheikhan for very 
valuable discussions. 
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